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ABSTRACT 

We have derived the temporal power spectra of the horizontal velocity of 
the solar photosphere. The data sets for 14 quiet regions observed with the G- 
band filter of Hinode/SOT are analyzed to measure the temporal fluctuation of 
the horizontal velocity by using the local correlation tracking (LCT) method. 
Among the high resolution (~0'.'2) and seeing- free data sets of Hinode/ SOT, we 
selected the observations whose duration is longer than 70 minutes and cadence 
is about 30 s. The so-called k-u diagrams of the photospheric horizontal velocity 
are derived for the first time to investigate the temporal evolution of convection. 
The power spectra derived from k-u diagrams typically have a double power law 
shape bent over at a frequency of 4.7 mHz. The power law index in the high 
frequency range is -2.4 while the power law index in the low frequency range is 
-0.6. The root mean square of the horizontal speed is about 1.1 km s _1 when we 
use a tracer size of C'4 in LCT method. Autocorrelation functions of intensity 
fluctuation, horizontal velocity, and its spatial derivatives are also derived in 
order to measure the correlation time of the stochastic photospheric motion. 
Since one of possible energy sources of the coronal heating is the photospheric 
convection, the power spectra derived in the present study will be of high value 
to quantitatively justify various coronal heating models. 

Subject headings: Sun: granulation — Sun: oscillations — Sun: photosphere 
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1. INTRODUCTION 

At least three types of convection are known to ubiquitously exist at the solar surface: 
granules, mesogranules, and supergranules. Their properties such as size, lifetimes, shapes 
have been studied by many authors (see reviews by, e.g., Leighton 1963, Spruit et al. 1990, 
Nordlund et al. 2009 and references therein). Through the magnetic field, the convective 
energy is transported upward to supply the sufficient energy to maintain the 1 MK corona. 
Therefore, the dynamics of the convective motions plays a key role in the so-called coronal 
heating problem, one of the most important issues in the solar physics. However, the 
temporal evolution of convection motion have not been sufficiently elucidated. 

The vertical motion in convection is easily observed by using Doppler analysis. This 
motion generates compressible waves such as magnetohydrodynamic slow mode waves or fast 
mode waves. Although a lot of energy can be transported upward by the waves caused by 
convection, compressible waves are not considered to contribute to coronal heating, because 



the wave energy will be significant 



or reflection in the chromosphere (IHollweg 



y reduced be: 



ore reaching the corona by shock dissipation 



19811 ). On the other hand, the horizontal motion 



in convection plays an important role in the coronal heati ng. The interaction betw een 



magnetic field and the convection generates Alfven waves (jUchida fc Kaburaki 



1974) while 



the stochastic phot ospheric moti on braids the field lines to store the energy in the corona 



as electric current ((Parker 



19831 ). Due to the lack of observations, it has been poorly 



understood what mechanisms contribute to coronal heating. 



The local correlatio n tracking (LCT) method 



horizontal velocity field ((November fe Simon 



1988 



commonly is used to derive the 



Berger et al. 



19981 ). Since LCT method 



uses apparent motion of granules to derive the velocity, it is better to use high spatial 
resolution and seeing-free data sets. From the horizontal velocity derived by LCT method, 
mesogranules and supergranules can be observed (e.g. Kitai et al. 1997, Ueno & Kitai 
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998, Shine et a 



Tit 



e et al. 



1989 



2000). Although various studies reveal the freq uency distributi on 



Berger et al 



19981 ) and the spatial power spectra (IRieutord et al. 



2000 



20081 ) of horizontal ve 



spectra ( jTarbell et al. 



ocity fields, few studies explicitly mention the temporal power 



19901 ) . In a previous study, we derived the photospheric horizontal 
velocity in 14 different quiet regions to justify the Alfven wave model for coronal heating 



( jMatsumoto fc Shibata 



20101 ) . In this study, we will show the temporal evolution of the 
horizontal motion of the photosphere in more detail. When discerning between various 
coronal heating models (see reviews by, e.g., Mandrini et al. 2000, Aschwanden et al. 2001, 
Klimchuk 2006, and references therein), the power spectra derived in the present study will 
be of high value. 



2. OBSERVATION AND DATA REDUCTION 

Hinode was launched on September 22, 2006 in order to investigate the unsolved 



problems of so 


ar p 


rysics, such as coronal heating. Solar Optical Telescope 


SOT) 


(Tsuneta et al. 


2008 


Suematsu et al. 


2008; 


Ichimoto et al. 


2008; 


Shimizu et al. 


2008 


), one 



of the three telescopes mounted on Hinode, is an optical telescope whose spatial resolution 
is unprecedentedly high (0 / /2-0 / /3, or 150-200 km) for a solar telescope in space. As opposed 
to ground based observations that are affected by atmospheric seeing, seeing-free data sets 
over a long time span can be continuously obtained by SOT. We selected 14 data sets of 
continuously observed quiet regions with the G-band filter between October 31, 2006 and 
December 29, 2007. Our basic selection criteria was for the data set to have a duration 
longer than 70 minutes, with a mean time cadence less than 32 s. Duration of the selected 
data sets ranges from 75 min to 345 min and the cadence is almost constant with the 
accuracy of 1 second. Our data sets are selected to be located within ± 100 arcsec from 
the disk center so that the line of sight is within 6 degrees from the local normal. Pixel 
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resolution of the CCD is 0'.'054 (40 km) and FOV is larger than 27" x 27"(20,000 km x 
20,000 km). 

We applied dark current subtraction and flat-fielding in the standard manner for all of 
the G-band images. The solar rotation, tracking error, and satellite jitter causes overall 
displacements between two images. Using cross correlation between two consecutive images, 
we have derived the optimum shift to adjust the overall displacements in an accuracy of sub 
pixel scale. 

Besides granular motion, p mode waves generally disturb the intensity in the G-band 
movies. In order to reduce t he contribution f rom p modes, we have applied so called 



sub-sonic filtering technique (ITitle et al. 



19891 ) to our data set. The sub-sonic filtering 



technique removes the power where u/k > V p h in Fourier space. In the present study, we 
have selected V P h = 5 km s -1 . 

The photospheric horizontal velocity of granular motion can be obtained by applying 
the LCT tec h nique to the G-band time series. In the LCT techniques developed by 



Berger et al. 



(119981 ). rectangular subfields or "tiles" are used, and the displacement of 
tiles between two consecutive images determines the horizontal velocity. Although the 
amplitude of the LCT velocity usually depends on the size of the tracer, we h ave fixed the 



tile siz e to be 0'.'4 in the present study, the same resolution as in the study of 



Berger et al 



(119981 ). From the horizontal velocity fields, divergence (V • V — d x V x + d y V y ) and rotation 
(V x V\ z = d x V y — d y V x ) can be derived, where x, y are the coordinates in the photospheric 
plane and z is the vertical coordinate. Figure [T] shows the G-band intensity image with 
horizontal velocity arrows (left), divergence (middle), and rotation (right), averaged over 
10 3 s. 
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3. 



RESULTS AND DISCUSSIONS 



We derived the apparent velocity in G-band data from Hinode/ SOT by using the LCT 
method. Figure |2] shows the so-called k-u diagram of the horizontal velocity normalized 
by the maximum power. The time series of LCT flow maps were Fourier transformed to 
estimate the power spectra. The dotted line in the figure represents u/k = 5 km s" 1 . Since 
we use the subsonic filtering method, spectral power above the dotted line is significantly 
reduced. The region in the lower frequency and the smaller wave number has greater power. 

Integrating k-u diagram over wave number space, we can estimate the frequency power 
spectrum density, P u . The power spectrum density is defined as, 



where V± is the LCT velocity, v min is the lowest frequency determined by total duration, 
and Umax is the highest frequency coming from the observational sampling time. The 
symbol, <>, denotes the temporal average. Figure |3] shows an example of the power 
spectrum density. Power spectra can be fitted to double power-law function, P v oc v ah when 
v < u\, (so-called break frequency) and P v oc v aa when v > u^. Typically, y/<V±>, z/ b , ol, 
and a.H of LCT velocity are 1.1 km s _1 , 4.7 mHz, -0.6, and -2.4, respectively in this work. 
By using the same analysis, temporal power spectra of divergence and rotation of velocity 
fields are also estimated (the dashed line and the dotted line in figure E]). The total power, 
ocl, and an of the divergence is 8.4 x 10~ 3 s^ 1 , -0.2, -1.5 respectively and the divergence 
will contribute to produce the acoustic mode waves. The total power, a^, and an of 
the rotation is 6.8 x 10~ 3 s _1 , -0.2, -1.3 respectively and the rotation is related to the 
generation of torsional Alfven waves. The power spectra of the LCT velocity has a harder 
power law index than that of its derivatives (V ■ V, V x V). 

The resulting power spectra are quite essential for evaluating models of coronal heating 




(1) 
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by Alfven waves ( 


Hollwcg 


1982 


Matsumoto & Shibat J 


201ok . L 



In their numerical sim ulations, nonlinear Alfv e n wav es are 



Matsumoto & Shibatal (120101 ) include 



driven by photospheric convection to heat the corona, 
the observed temporal power spectra shown here to generate Alfven waves and succeed in 
producing sufficient energy flux to heat the corona. 

Next, in order to derive the "correlation time" of the stochastic convection motion, we 
derived the autocorrelation function (ACF), 

</(*)/(* + r)> 



R(-, 



< f(t) 2 > 

, where / is one of physical variables associated with turbulent motion and <> denotes 
temporal averaging process. The correlation time is determined as the e-folding time 
of the ACF. Correlation time often scales with coronal heating rate in various heating 



(2) 



models that use stochastic c onvection motion (jSturrock fc Uchida 



1981 



Parker 



Galsgaard fc Nordlund 



19881 : 



19961 ). Figure H] shows the ACF of the G-band intensity fluctuation 
(solid line), horizontal velocity (dashed line), divergence (dotted line), and rotation 
(dash-dotted line). For G-ba nd intensity the co rrelation time is about 250 sec, which is 



comparable with the value in 



Title et al. 



(119891 ) and references therein. The horizontal 



velocity has sh o rter c orrelation time (~ 100 sec) than the G-band intensity fluctuation. 



Rieutord et al. 



(120001 ) also derived the ACF of the horizontal velocity and found a 
relatively larger correlation time of about 35-45 min, because they were looking at larger 
spatial/temporal scale phenomena. ACFs of rotation and divergence fall rapidly and their 
correlation times are about 50 sec. The correlation time in the coronal heating models is 
usually defined by the velocity correlation time (e.g. Sturrock & Uchida 1981). However, 
the lifetime of granules, which corresponds to the intensity correlation time, is usually used 
as the correlation time in the previous studies. Since the velocity correlation time turned 
out to be shorter than the intensity correlation time, the coronal heating rate will be greatly 
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reduced. 



Let us explain why the correlation time of intensity fluctuation (rj) is longer than that 
of velocity (ry) qualitatively. The power spectrum of intensity (Pj) can be regarded as the 
power spectrum of spatial scale. From the dimensional analysis, the power spectrum of 
velocity (Py) is proportional to Pju 2 so that Pj has a spectrum that is at most a power 2 
softer than Py. Since Fourier transform of ACF corresponds to power spectrum, softening of 
the power spectrum means hardening of the ACF decay or long correlation time. Therefore 
it is reasonable that tj is longer than ry. From the similar dimensional analysis, the shorter 
correlation time of divergence or rotation can be explained. 

If we integrate k-u diagram over frequency space, the spatial power spectrum density 
can be estimated (Figure E]). The power law index for the high k range (k ~ 10 -3 ) becomes 
-5/3, which repr esents the Kolmqgoroy type turbulence. Our results are complementary 



to the results of 



Espagnet et al. 



Rieutord et al 



3) 



(120081 ) since they investigate large spatial structures. 



(119931 ) showed the same power law index in the power spectra of vertical 



photospheric velocity. These results indicate that the granules are Kolmogorov turbulent 
eddies. 

We also analyze motion of test particles in the LCT flow field. From the spatial 
deviation of test particles (5) and the elapsed time (r), diffusion coefficient (D ph = 6 2 /2r) 
can be derived, if we assume motion of each test particle can be approximated by random 



walk. For our data set, 
comparable with that of 



;he diffusion coefficient D p h is around 500 km 2 s \ which is 



Tarbell et al. 



current cascading coronal heating model (IVan Ballegooijen 



(1990 ). When we apply t 



lis value to the electric 



19861 ) . heating rate that scales 



with the diffusion coefficient of the photospheric horizontal motion becomes less than half 
of the required rate. Therefore when we consider coronal heating in the quiet sun, current 
cascading may not contribute significantly to the heating rate. 
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Since the LCT method tracks the apparent motion of the photospheric convection, the 
velocity derived here would contain some artifacts from the apparent velocity even though 
we have carefully removed the effect of acoustic waves by subsonic filtering method. The 
decrease in density or increase in tem perature causes the destruction of CH molecules, 



resulting in G-band intensity changes (jSteiner et al. 



200 ll ). These intensity changes can 



affect the LCT velocity. It is important to compare the result of LCT measurements 



with the recent realistic 3D s imulations (jStein fc Nordlund 



Georgobiani et al 



2006 



20071 ). 



1998 



Rieutord et al 



2001 



In conclusion, k-u diagrams of LCT velocity are derived for the first time to investigate 
the temporal evolution of the photospheric convection. Integrating the power spectra 
over wave number space generally reveals a double power law spectral shape whose break 
frequency is about 4.7 mHz. The power law index in the low frequency range is -0.6 while 
the power law index in the high frequency range is -2.4. The correlation time of G-band 
intensity and horizontal velocity is about 200 sec and 100 sec respectively while that of 
divergence and rotation is 50 sec. The diffusion coefficient of photospheric convection is 
around 500 km 2 s _1 , which is not sufficient to heat the corona above quiet region through 
current cascading. By using the power spectra derived in the present study, some of the 
coronal heating models will be justified (e.g. Matsumoto & Shibata 2010). 



Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
domestic partner and NASA and STFC (UK) as international partners. It is operated by 
these agencies in co-operation with ESA and NSC (Norway). Scientific operation of the 
Hinode mission is conducted by the Hinode science team organized at ISAS/JAXA. This 
team mainly consists of scientists from institutes in the partner countries. Support for the 
post-launch operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, 
ESA, and NSC (Norway). Takuma Matsumoto is supported by the Research Fellowship 
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Fig. 1.— G-band image with horizontal velocity arrows (left), divergence (middle), and 
rotation (right). The velocity field and its derivatives are averaged over 10 3 s. 
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l io inn 
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Fig. 2. — k-u diagram of horizontal velocity fluctuation. The dotted line indicates u> = 5 
km s _1 k. The gray scale color represents the power spectrum density normalized by the 
maximum power in log scale. 




Fig. 3. — Temporal power spectrum density of the photospheric horizontal velocity (solid 
line), rotation (dotted line), and divergence (dashed line). After the horizontal velocity is 
derived by the LCT technique, the power spectrum is estimated by FFT. 
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Fig. 4. — Normalized autocorrelation function of G-band intensity fluctuation (solid line), 
horizontal velocity (dotted line), rotation (dashed line), and divergence (dash dotted line) 
defined by the equation [2j 
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Fig. 5. — Spatial power spectrum density of the photospheric horizontal velocity. The power 
law index of the dotted line is -5/3. 
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